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Abstract 

Density functional study of strain effects on the electronic band structure and transport prop- 
erties of the graphene nanoribbons (GNR) is presented. We apply a uniaxial strain (e) in the x 
(nearest-neighbor) and y (second nearest-neighbor) directions, related to the deformation of zigzag 
and armchair edge GNRs (AGNR and ZGNR), respectively. We calculate the quantum conduc- 
tance and band structures of the GNR using the Wannier function in a strain range from —8% to 
+8% (minus and plus signs show compression and tensile strain). As strain increases, depending 
on the AGNR family type, the electrical conductivity changes from an insulator to a conductor. 
This is accompanied by a variation in the electron and hole effective masses. The compression e x 
in ZGNR shifts some bands to below the Fermi level (Ef) and the quantum conductance does not 
change, but the tensile e x causes an increase in the quantum conductance to 10e 2 //i near the Ef. 
For transverse direction, it is very sensitive to strain and the tensile e y causes an increase in the 
conductance while the compressive e y decreases the conductance at first but increases later. 

PACS numbers: 71.70.Fk, 72.80.Rj, 73.61.Wp 
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I. INTRODUCTION 



The discovery of graphene in 2004 and its potential application in nanotechnology has 
attracted considerable interest 1 3 . The unique properties of graphene due to massless Dirac 
fermions, high surface to volume ratio, high crystal quality in 2D, and high mechanical 
strength recommend it as a promising material for technologies such as nanoelectromechan- 
ical systems (NEMS) and nanosensor^H3 < u se f graphene in nanodevices is accompanied by 
induced strains and stresses. These may originate from phonon-induced lattice vibrations, 
lattice mismatched film growth, or applied external stress, and it is therefore necessary to 
take them into account. 

Recently, research has focused on strained graphene in order to study the effect of strain on 
the electronic structure propertied^^R Deformations in a graphene sheet affect its electron 
transport by decreasing its conductance at low densities, making electrons behave in the 
same way as in a nanoribbon or a quantum dot 10 . Theoretical studies of the strain effects, 
followed by recent experiments, have attempted to understand and control the interplay 
between the strain and the electrical transport in graphenJ^. 

In a nanoribbon, the wave function is localized in two dimensions. Because of the finite 
size of a GNR, Bloch orbit als cannot be used directly in electronic transport calculations, 
and quantum conductance should be computed via lattice Green's function and localized 
orbital representation of the electronic states in real space. Here, we employ maximally 
localized Wannier functions (MLWF) to compute the electrical conduction propertied. 

In this paper we perform an ab initio study of the strain effect on a flat GNR as the 
structural motif for carbon nanotubes, and other graphene-based material^. We calculate 
the coherent transport properties of nanostructures from first principles by combining them 
with the Landauer approacbPH The connection between first principles and Landauer ap- 
proach is provided by using MLWF representation. Using MLWFs, we calculate the band 
structure and the quantum conductance of a GNR under various strains. The results show 
that the transport properties of a GNR depend on the strain. As the strain increases, de- 
pending on the AGNR family type, the electrical conductivity changes from an insulator to 
a conductor, and this is accompanied by a variation in the electron and hole effective masses. 
The compression e x in the ZGNR shifts some bands to below the Ef and the quantum con- 
ductance does not change, but the tensile e x causes an increase in the quantum conductance 
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to 10e 2 /h near the Ef. For transverse direction, the ZGNR is very sensitive to the strain, 
and the tensile e y causes an increase in the conductance while the compressive e y decreases 
it at first, but increases later. 

II. COMPUTATIONAL DETAILS 

First-principles DFT-based calculations, which have been successfully applied to the 
study of graphen^O, have been performed in this paper. All our DFT calculations were car- 
ried out with the PWSCF cod d 17 | 1S | , within the local density approximation (LDA) and using 
the ultrasoft Rabe-Rappe-Kaxiras-Joannopoulos pseudopotentials 19 . We use an MLWF ba- 
sis set which can be obtained by a unitary transformation of the occupied ground-state plane 
wave eigenfunctions of the density functional theory with the cutoff energy of 340 eV. The 
total system consists of 12 and 18 C atoms (graphene supercell) in a slab geometry in two 
dimensions (Fig. [I]) and with a distance of 14 A between the adjacent layers. The sampling 
of the Brillouin zone is done using a 1 x 1 x 8 Monkhorst-Pack gricP^, which is tested to give 
converged results for all properties that are calculated. MLWFs are computed following the 
method of Marzari and Vanderbilt (MV) using wannier90 cod d 21 1 221 . For entangled energy 
bands, the method of Souza, Marzari and Vanderbilt (SMV) is usecpl. The coordinate op- 
timization is done via conjugate gradient molecular dynamics with geometry constrains in 
induced strain direction, and the maximum atomic displacement is limited to 0.1 A for each 
molecular dynamics step in order to preserve the stability of calculations. 

III. RESULTS AND DISCUSSION 

The structure studied in this work is made from a graphene sheet that is sheared along 
a zigzag or an armchair edges and passivated by hydrogen atoms (Fig. [I]). Recently, the 
effects of edge and size on the GNR band structures and also the effect of charge transfer 
from substrate to a graphene layer have been studiecpSEH. Here, we have neglected the 
charge transfer from substrate, and the edge and size effects, and have only considered the 
strain effects. We represent the deformation of the GNR by the strain defined as Aa^/a^ 
where is the strain induced superlattice constant in the i — x,y directions. To study the 
strain effects, we consider the longitudinal and transverse strains separately i.e., in the x 
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(a) 



(b) 



FIG. 1: Graphene nanoribbon supercell. 

and y directions. 

First we calculate the electronic structure of the AGNR at zero strain. An AGNR with 
widths w — 9, 10, and 11 is chosen to represent three typical families corresponding to 3n, 
3n+l, and 3n+2 (Fig. [I]), respectively, similar to the previous work^SESl We present every 
family with w key as w-AGNR. Fig. [2] presents the band structures, the density of states 
and the quantum conductance of a 9- AGNR. Neglecting the edge effects, due to passivation, 
Van der Waals bonds, that create band dispersion in out of the basal plane, are responsible 
for band mixing. The band structure is spanned by the bonding combination of sp 2 orbitals 
(the graphitic backbone) and p z orbitals (jr and 7r* bands). The hopping integral of the 
7r bands dispersion near the Ef determines the band gap, the transport property, and the 
effective mass of the charge carrier. As shown in Fig. [2j at T point, there is a direct gap of 
1.4 eV between the valence and conduction bands. In the valence band, the heavy hole band 
(Vi subband) is separated from the light hole band (V2 subband). Van Hove singularity in 
the density of states shows this band splitting. 

Using retarded and advanced Green's functions, based on the Landauer formalism, quan- 
tum conductance is computed. MLWFs, evaluated by a combination of plane- wave electronic 
structure, is applied for further conductance calculation. To construct the wave functions 
we select an energy window with E G [—6, 6] eV around the Ef (taken as the reference zero). 
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This energy window contains all the occupied and empty states. Fig. [2] shows the quantum 
transport of the 9-AGNR and its value near the Ef is zero. Here, conductance includes sum 
over all transmission possibilities that are between and 1. Therefore, there are no channels 
to connect the two sides of the 9-AGNR near the Ef. 

Fig. [3] shows the band structures of three families of AGNR under various uniaxial tensile 
(s > 0) and compressive (e < 0) strains. The band structures of the AGNR families are 
sensitive to the strain, and the band gap changes non-symmetrically for the tensile and the 
compressive strains as shown in FigjiJ For the 3n family in Fig, 4(a) the tensile e y and the 
compressive e x increase the band gap at first and then reach a maximum value at a strain of 
about e x = —0.02 and e y = 0.02, and finally the application of further strain decreases the 
band gap. Band gap variation for the 3n+l group is shown in Fig. |4(b)[ It suppressed to a 
narrow gap when the induced strain reaches e y = 0.04 and e x = —0.02, while the compressive 
strain in y direction and the tensile strain in x direction increase the band gap to above 1 
eV. Band gap variation for the 3n+2 family is completely different, and has a maximum 



and a minimum (Fig. 4(c)). The tensile e y decreases the gap, while the compressive strain 
causes it to increase the gap. The corresponding maximum gap occurs at e x = 0.02 and 
e y = —0.02 and the minimum is at e x = —0.04 and e y = 0.04. As shown in Fig. [4j all the 
behavior in the y direction is inverse of that in the x direction, and our results demonstrate 
that the band gap does not vanish under any condition. 

Another result from Fig. 3(a)| and |3(b) is the variation in the effective mass of fermions 
as the strain induced curvature of bands at Y point, that is related to the effective mass, 
changes. When the tensile strain becomes compressive, in addition to variation in the 
effective mass, the heavy hole band near the band edge (Vi subband) exchanges with the 
light hole band (V2 subband). The variation in the effective mass is different for w =9 and 
10, 11. The deformation of the 9-AGNR only causes a variation in the effective mass, but for 
the 10 and 11- AGNR family it can decrease the fermion's effective mass to near zero (Fig. 
3] for s x = —0.04 and e y = 0.04). In this case, dispersion relation near the Ef is converted 
to quasi-linear by a narrow gap. In addition, band deformation, which depends on strain, 
results in two types of excitons; first type relates to heavy electron-hole, and the second 
corresponds to light electron-hole. 

Figj5] shows the quantum conductance of the 9, 10 and 11-AGNR under various induced 
e y and e x at T = K. Only electrons at Ef play a role in the quantum conductance. The 9 
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FIG. 2: Band structures, density of states and quantum conductance of a 9-AGNR at zero strain. 

and 10-AGNR, at every induced strain, have zero conductance but for the 11-AGNR there 
are two special strains {e x = -0.04 and e y = 0.04) wherein the ribbon conductance is nonzero 
due to a narrow gap. At finite temperature, every subbands near the Ef can contribute to 
the quantum conductance. Current for a finite bias voltage (V) is given 



1 = [ d i [f{£ + "> - m]G{£) * r d > d i ]a{£) = v r d£ fe^ 



e) 



where / is the Fermi-Dirac distribution function and \± is the chemical potential. At a finite 
temperature /i |^ has a Gaussian peak that at T = K becomes a Dirac delta function. 
We estimate that at room temperature the current can be approximated by integrating over 
£ G [—0.3, 0.3]. Fig. [6] shows the conductance (G = ^) of an AGNR family versus strain at 
room temperature. Because of relatively high band gap of the 9-AGNR, it has a relatively 



low conductance at room temperature for a wide range of strains (Fig 6(a)). However, the 



10 and 11-AGNR have transition between two conductance levels for the compressive and 



the tensile strains. Fig. 6(b) shows the conductance of a 10-AGNR that was computed from 
Eq. [TJ For e x > and e y < 0, the AGNR is an insulator while for other strain values it is 



a conductor. As shown in Fig. 6(c) the conductance of the 11-AGNR for strain values of 
s y = —0.02 and e x = 0.02 decreases due to the increase of the band gap. 

We now present the electronic structure of an ZGNR at zero strain. Fig. [7] presents the 
band structures, the density of states and the quantum conductance of an ZGNR. The band 
structure near the Ef has a linear dispersion and the group velocity near this level is about 
2.1 x 10 6 ms~ l . To calculate the conductance and the density of states we construct the 
wave functions in the energy window with E G [—6, 6] eV around the Ef that is taken to be 




(a): Strain in x direction for three families of 
ribbons 




(b): Strain in y direction for three families of 
ribbons 



FIG. 3: Band structures of n-AGNR for three families of ribbon under various strains. 
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FIG. 4: Energy gap of an AGNR under various strains for w = 3n, 3n+l and 3n+2. 

the reference zero. Quantum transport of the ZGNR near the valence band edge is about 
6e 2 /h. 

Fig. [8] shows the band structure of a ZGNR under various uniaxial strains. The induced 
e x cause an increase in the group velocity of the Dirac fermions from 2.0 x 10 6 ms~ l to 



maximum 2.3 x 10 6 ms~ l for the strain from -0.08% to 0.08% (Fig. 8(a)), but the induced 
£ y cause an decrease in the group velocity from 2.7 x 10 6 ms~ l to minimum 2.0 x 10 6 ms~ l 
for the same strain range (Fig. |8(b)[ ). The variation in the band structures influence the 



electrical conductance of the GNR. Fig. |8(c) and |8(d)| show the quantum conductance of 



the ZGNR under various e x and e y . By comparing these results with those related to the 
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(a): The Landauer quantum conductance variation (b): The Landauer quantum conductance 
in the x and y direction. variation in the x and y direction. 

X 11-AGNR Y 



5^ -0.08 

L - J — 



£ 5 
^0 

(N 

o 5 


5 




04 



0.04 




-0.08 



0.08 



L j 



-0.04 



0.04 



0.08 



2-1012-2-1012 
E-E f (eV) E-E f (eV) 

(c): The Landauer quantum conductance variation 
in the x and y direction. 

FIG. 5: Energy gap of an AGNR under various strains for w = 3n, 3n+l and 3n+2. 



relaxed ZGNR (Fig. [7]) we find that the strain can increase the conductance of the ZGNR 
up to two times. The compressive e x shifts some bands to below the Ef but doesn't change 
the quantum conduction near the Ef, while the tensile e x enhances it to 10e 2 //i. In the y 
direction, the ZGNR is very sensitive to the strain and the tensile e y causes an increase in 
the conductance, while the compressive e y decreases it at first but increases later. 
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FIG. 6: Conductance of an AGNR under various strains for w = 3n, 3n+l and 3n+2 at room 
temperature. 
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FIG. 7: The band structures, the density of states and the quantum conductance of an ZGNR at 
zero strain. 



IV. CONCLUSIONS 



Density functional studies of strain effects on the electrical transport properties of the 
graphene nanoribbons are presented. By applying a uniaxial tensile strain in the x and y 
directions, the electronic properties of the graphene nanoribbons were studied. Using the 
Wannier functions, the band structure and density of states were calculated for different 
strains from —8% to 8%. It is observed that as the strain increases, depending on AGNR 
family type, the electrical conductivity changes from an insulator to a conductor, and this 
is accompanied by the variation in the electron and hole effective masses. Compressive e x 
in the ZGNR shifts some bands to below the Ef and the quantum conductance does not 
change but the tensile e x causes an increase in the quantum conductance to 10e 2 //i near 
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(c): The Landauer quantum conductance (d): The Landauer quantum conductance 
variation by x direction strain. variation by y direction strain. 



FIG. 8: The band structures and the quantum conductance of an ZGNR under various strains. 

the Ef. In the transverse direction, the ZGNR is very sensitive to the strain and the tensile 
£ y causes an increase in the conductance while the compressive e y decreases it at first but 
increases later. 

Acknowledgments 

H. R-T. acknowledges the support from Iran National Science Foundation. We thank Dr. 
M. Farjam for useful discussions. 



11 



* Electronic address: rani-tabar@n ano.ipm.ac.irl 

1 A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 
81, 109 (2009). 

2 A. Reina, X. Jia, J. Ho, D. Nezich, H. Son, V. Bulovic, M. Dresselhaus, and J. Kong, Nano. 
Lett. 9, 30 (2009). 

3 K. Kim, Y. Zhao, H. Jang, S. Lee, J. Kim, K. Kim, J. Ahn, P. Kim, J. Choi, and B. Hong, 
Nature 457, 706 (2009). 

4 A. K. Geim and K. S. Novoselov, Nat. Mater. 6, 183 (2007). 

5 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grig- 
orieva, and A. A. Firsov, Science 306, 666 (2004). 

6 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, and 
A. Dubonos, S. V .and Firsov, Nature 438, 197 (2005). 

7 K. S. Novoselov, E. McCann, S. V. Morozov, V. I. FalKo, M. I. Katsnelson, U. Zeitler, D. Jiang, 
F. Schedin, and A. K. Geim, Nat. Phys. 2, 177 (2006). 

8 E. Cadelano, P. L. Palla, S. Giordano, and L. Colombo, Phys. Rev. Lett. 102, 235502 (2009). 

9 V. M. Pereira, A. H. Castro Neto, and N. M. R. Peres, Phys. Rev. B 80, 045401 (2009). 

V. M. Pereira and A. H. Castro Neto, Phys. Rev. Lett. 103, 046801 (2009). 

1 M. Farjam and H. Rafh-Tabar, Phys. Rev. B 80, 167401 (2009). 

2 W. Bao, F. Miao, Z. Chen, H. Zhang, W. Jang, C. Dames, and C. Lau, Nat. Nanotechnol. 4, 
562 (2009). 

3 I. Souza, N. Marzari, and D. Vanderbilt, Phys. Rev. B 65, 035109 (2001). 

4 R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical properties of carbon nanotubes (Im- 
perial College Press, 1998). 

5 M. B. Nardelli, J. L. Fattebert, and J. Bernholc, Phys. Rev. B 64, 245423 (2001). 

6 O. Leenaerts, B. Partoens, and F. M. Peeters, Phys. Rev. B 77, 125416 (2008). 

7 S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Giannozzi, PWscf (Plane-Wave Self- Consistent 
Field) codes. 

8 P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. Chiarotti, 
M. Cococcioni, I. Dabo, et al., J. Phys. Condens. Mat. 21, 395502 (2009). 



12 



A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. Joannopoulos, Phys. Rev. B 41, 1227 (1990). 
H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 
N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 (1997). 

A. Mostofi, J. Yates, Y. Lee, I. Souza, D. Vanderbilt, and N. Marzari, Comput. Phys. Commun. 
178, 685 (2008). 

F. Varchon, R. Feng, J. Hass, X. Li, B. N. Nguyen, C. Naud, P. Mallet, J. Y. Veuillen, C. Berger, 
E. H. Conrad, et al., Phys. Rev. Lett. 99, 126805 (2007). 

D. Finkenstadt, G. Pennington, and M. J. Mehl, Phys. Rev. B 76, 121405 (2007). 

Y. W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 97, 216803 (2006). 

Z. F. Wang, Q. Li, H. Zheng, H. Ren, H. Su, Q. W. Shi, and J. Chen, Phys. Rev. B 75, 113406 

(2007). 

V. Mujica, M. Kemp, and M. Ratner, J. Chem. Phys. 101, 6849 (1994). 



13 



